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Abstract 

We demonstrate the level statistics in the vicinity of the Anderson transi- 
tion in d > 2 dimensions to be universal and drastically different from both 
Wigner-Dyson in the metallic regime and Poisson in the insulator regime. The 
variance of the number of levels N in a given energy interval with (N) S> 1 
is proved to behave as (iV) 7 where 7 = 1 — (vd)~ 1 and v is the correlation 
length exponent. The inequality 7 < 1, shown to be required by an exact 
sum rule, results from nontrivial cancellations (due to the causality and scal- 
ing requirements) in calculating the two-level correlation function. 
PACS numbers: 71.30. +h, 05.60,+w, 72.15.Rn 

Typeset using REVTeX 



*To be published in PRL, vol. 72, p.888 (1994) 



1 



The problem of level statistics in random quantum systems is attracting considerable 
interest even now, four decades after the pioneer works of Wigner and Dyson This is 
because of the universality of the Wigner-Dyson statistics which makes it relevant for a large 
variety of quantum systems 0. 

For the problem of a quantum particle in a random potential, the Wigner-Dyson statistics 
is known to be applicable for finite systems in the region of extended states which will be 
referred to as a metallic region. With increasing the random potential, the system undergoes 
the Anderson transition into the insulator phase ||, where all states are localized. In this 
region, the statistics of energy levels is expected to be Poisson. 

There is, however, the third region, namely, the critical region in the vicinity of the An- 
derson transition where the spectral statistics is believed to be still universal |7]|| , although 
different from both Wigner-Dyson and Poisson. As the critical region can not be consid- 
ered perturbatively or semiclassically, nearly nothing is known about the third universal 
statistics. 

The first attack at this problem has been done in Ref. |7|] where the simplest statistical 
quantity, the variance ((5N) 2 ) of the number of energy levels N in a given energy interval of 
the width E, has been considered (SN = N — (N), and (. . .) denotes the ensemble average 
over the realizations of the random potential). The dimensional estimation made in Ref. |7j 
has resulted in ((5N) 2 ) = a(N), thus being different from the Poisson statistics only by a 
certain number a < 1. 

We will see, however, that this result contradicts to an exact sum rule resulting from the 
conservation of the total number of levels. The point is that in the dimensional estimations |7j 
analytical properties of diffusion propagators have not been taken into account. We will show 
that the analytical properties resulting from causality together with certain scaling relations 
near the Anderson transition make the (iV)-proportional contribution to the variance to 
vanish. We will calculate the spectral density correlation function and deduce from it the 
following universal relationship between the variance and the average number of levels in 
the energy interval E 
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that holds exactly at the mobility edge. Here v is the correlation length exponent, and the 
factor ad/3 is universal in a sense that it is determined completely by the dimensionality d 
and the symmetry class of the Dyson ensemble {(3 = 1, 2, or 4 for unitary, orthogonal, and 
symplectic ensembles, respectively). For many systems v « 1 so that 7 ~ 2/3 for d = 3. In 
general, Eq. (1) suggests a new way of determining v. 

This is the main result of the paper. It demonstrates that in the vicinity of the Anderson 
transition there really exists the third universal statistics. It governs the spectral fluctuations 
that are weaker than for the Poisson statistics, ((<5iV) 2 )p ~ (N), but much stronger than for 
the Wigner-Dyson statistics, ((<5iV) 2 ) w _ D ~ In (AT). 

All the three statistics are universal and exact in the same limit: 

L -> 00 E/A = (N) = const > 1 (2) 

where L is the sample size. In this limit, the mean level spacing A = (uoL d )~ 1 tends to zero 
(uq is the mean density of states), but the number of levels in an interval E is kept finite, 
although very large. 

The new level statistics describes the fluctuations in an energy band \e — Eq\ < E/2 
centered exactly at the mobility edge Eq = e c . For the critical regime to be achieved the 
correlation length L c (e) which diverges as \e/e c — l\~ u must exceed the sample size for all e 
in the energy band E. Due to this uncertainty L c = (E/e c )~ p ', and 

LJL = (N)~»{L/\) d »-\ (3) 

where A = (z/o£ c ) _1//d - Then the Harris criterion v > 2/d ensures L c /L — ► 00 in the 
limit (0) for the energy band centered at e = e c- in the same limit, the Wigner-Dyson and 
Poisson statistics describe exactly the fluctuations in bands centered at Eq > e c (the metallic 
region) and Eo < e c (the insulating region), respectively. The limit (0) is required, therefore, 
to avoid mixing the levels belonging to different regions as well as to make the finite-size 
corrections vanishing. 
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We consider the spectral density correlation function 

R(u) = —(v(e)v{e + u))-l, (4) 

where u(e) is the exact density of states at the energy e. Note that the function R(u) has 
a singular term S(u) resulting from the self-correlation of energy levels. 

Before deriving the announced result, Eq. ([!]), we demonstrate that the exact sum 
rule prohibits the variance ((5N) 2 ) to be (iV)-proportional. The conservation of the to- 
tal number of energy levels for any non-singular random potential may be written down as 
f^oo[ v { £ + u ) ~ v o]duJ = 0. It leads to the sum rule: 

/oo 
R{s)ds = 0, s = u/A. (5) 
-oo 

The variance ((5N) 2 ) of the number of levels N in the energy band of the width E centered 
at a certain energy e (e.g., at the Fermi level e F ) is given by 

/ \ r( N ) 

((5N) 2 ) = l w ((N)-\s\)R(s)ds. (6) 

Then 



d((SN) 2 ) 
d(N) 




If the function R(s) is universal in a sense that it does not depend on any parameter, then 
the only condition (N) ^> 1 is sufficient, due to the sum rule (|5]), to make the integral in 
the r.h.s. of Eq. (|7|) to be arbitrary small. Therefore, in this case {(SN) 2 )/ (N) — > 0. 

The universality assumption is crucial for vanishing the contribution to the variance 
proportional to (N), or the higher power of (N). However, a finite disordered sample is 
characterized by a set of relevant energy scales that obey in the metallic limit the following 
inequalities: 

A < 1/td < l/r < e F (8) 

where td = L 2 /D is the time of diffusion through the sample, D is the electronic diffusion 
coefficient in the classical limit, D = v 2 r/d, r is the elastic scattering rate, h—1. Naturally, 
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for sufficiently large (N) the function R(s) depends not only on s. It results in ((6N) 2 ) oc 
(tdA) 3 / 2 (N) 3 ^ 2 in an energy band of the width E ^> 1/td ||. We will show elsewhere 
that higher than (N) contribution arises also in the critical region (A~1/td <C 1/t~£ f ) 
where it is proportional to (N) 1+a (rA) a (0 < a < 1 is a certain critical exponent). Both 
these nonuniversal contributions could be of importance for finite systems. However, they 
do vanish in the limit (|]). In this limit only the universal contributions to the variance 
survive. 

In the insulating regime, the above speculations are not applicable for estimating the 
integral in Eq. (0). The reason is the existence of the additional energy scale Ag = l/Vo£ d = 
A(L/£) d which is a typical spacing for states confined to a localization volume £ d centered 
at some point. Since such states are repelling in the same way as extended states in metal 
confined to the whole volume L d , the function R(s) at s ^ is expected to be similar 
to the Wigner-Dyson function Rw-d(w/ A) with substituting A by A^. Such a function 
R(s) = (£/ L) d R w _ £>[(£/ L) d s], which obviously obeys the sum rule @, is not universal at 
all scales and reduces to a constant — (£/L) d for s <C Therefore, in the limit (0) the 

regular part of R(s) makes no contribution to the r.h.s. of Eq. (|7|). Then d((5N) 2 ) /d(N) is 
exactly equal to 1 due to the singular <5(s)-term in R(s). 

Now we turn to microscopic calculations. In the metallic region, R(u) is given by the 
two-diffuson diagram U that is convenient to represent (see for detail Ref. 0]) as in Fig. 
1(a), separating the diffusion propagators (wavy lines). Both in the metallic region for 
cu^A and at the mobility edge one should consider also 2n-diffuson corrections (Fig. 2(a) 
for n=2). In all diagrams, the polygons with 2n + 1 vertices are made from the electron 
Green's functions that decrease exponentially over the distance of the mean free path. Thus, 
all vertices of any polygon correspond to the same spatial coordinate and its ensemble- 
averaged contribution reduces to a constant which we denote ^o r2 ™X2n+i, where X2n+i are 
dimensionless complex numbers. Then the general expression for the 2n-diffuson diagram 
in the momentum representation is given by 
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A2n| 12 ( 2n 1 

^»n= lx 2 2 : +l] ReE kll^q,) • (9) 

71 P <?l"<?2n[ j = l J 

Here P(u, q) is the exact diffusion propagator, 5 q = 5 gi+ ... +g2ni0 , and the factor accounts 
for the number of diagrams in different ensembles where some channels of propagation are 
suppressed. 

In the metallic region, P(u,q) = (Dq 2 — iuo)~ l . For ujt£, <C 1 (that corresponds to 
t ^> To), the excess particle density is distributed homogeneously over the whole sample 
so that only q = contribution of each diffuson survives in Eq. (|9|). For u ^> A, only the 
two-diffuson diagram (n — 1) is relevant || so that (with \3 — i) 

R( s ) = -(tt 2 /?* 2 )- 1 . (10) 

The sum rule (0) allows to calculate ((5N) 2 ) in the energy interval where (N) ^> 1 
using only the perturbative result fllOD. One represents the first term in Eq. (||) as 
—2 (JV) -R(s) (is which is a constant of order 1. The second term in Eq. (||) diverges 
only logarithmically. Restricting it to the perturbative region with a cutoff at s*sl, one 
reproduces the Wigner-Dyson result Jl| with the accuracy up to a constant of order 1: 

((6N) 2 } = -L ln(iV). (11) 

For E ^> 1/td, this result does not hold in the metal [g| where such a width is unreachable 
in the universality limit ([|), though. On the contrary, at the mobility edge A ~ 1/td and 
any interval with (iV) ^> 1 has the width E ^> 1/r^. That is why one expects the variance 
((5N) 2 ) to deviate drastically from that in Eq. (|TT|). 

At the mobility edge, P(u, q) may be expressed as 

P(w,q) = [D(cu,q)q 2 -iu]' 1 . (12) 

which is the most general expression compatible with the particle conservation law. Although 
the exact diffusion coefficient D(u, q) here is unknown, the scaling and analytical properties 
of the diffusion propagator enable us to determine R(s) for s ^> 1. Since the propagator 



P(t, r — r ; ), that is the space-time Fourier-transform of P(u,q), is nonzero only for t > 
(causality) and real, P(u, q) is analytical in the upper half-plane of the complex variable uj 
and satisfies the relation P*(uj,q) = P(—u,—q). Using also the spatial isotropy, one has 
P*(u,q) = P(-u,q). 

At the mobility edge in the limit L — > oo, the scaling arguments allow to express P(uj, q) 
in terms of the dimensionless scaling function F depending on qL^, the ratio of the only two 
lengths characterizing the system [11]]. Here L u is a characteristic length of the displacement 



of a diffusing particle for the time uj x . At the critical point, a dimensional estimation yields 

Luj ~ (ujv y l ' d . (13) 
With the standard definition L u = \D(uj) / uj\ l l 2 , Eq. ( |13|) reproduces the well-known scaling 
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D(o;) oc L 2 J d oc o; 1 " 2 ^. (14) 

Using the scaling relation ([IB]) , we obtain 

p( w , g ) = (-ia;)- 1 ^ (z) , 2 = -^z/ g" d , (15) 

where due to the above analyticity requirements, uj contains an infinitesimal imaginary part, 
and the function F(z) is analytical for Re^ > and satisfies the condition 

F*(z) = F(z*). (16) 



In the static limit P(u — > 0,q) oc q~ d at the critical point JT3|. In the opposite limit, 
L u q <C 1, the diffusion propagator has the form (|l|) with the diffusion coefficient (p~4[) 
depending only on u>. That results in the asymptotics 



(17) 




where are real coefficients of order 1. 
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Now we substitute Eq. (|TJ) into Eq. @, change S g .for L d J d d qj/ (2ii) d , and represent 
5 q as / d d r exp(ir X) q^) . Dividing the integration over r into that over the surface (Sj) and 
radius of the d- dimensional sphere, and introducing dimensionless variables kj = qjr and 
( = ujvtfr d , we reduce Eq. ([)]) to 

A(-iyWil 2 fjo * /-^ kjCOSd] 



-^ 2n ~~ a 2a / II / , o 

[5n A d J jJiJ (2ir) 



In 

/ dcr 2n Ren^(-<^)- ( 18 ) 

J0 J=l 



X 



A dimensional estimation of this integral would give A/a; so that R(s) ~ 1/s. Having sub- 
stituted this into Eq. (|), one would obtain ((SN) 2 ) ~ (N) In (N) which is strictly prohibited 
by the sum rule, as shown above. 

However, it follows from Eq. flT6f) that the real part of the product of F-functions in Eq. 
(plf) is an even function of (. Thus, the integration over £ can be extended to the whole real 
axis. Taking into account the asymptotics ([TTf ) and the analyticity of the function F(—iQ in 
the upper half-plane of the complex variable (, one concludes immediately that the integral 
(p"8|) equals zero. 

Therefore, R(u) = for u ^> A in the limit (0). For large but finite L, one has to consider 
corrections to the diffusion propagator proportional to powers of the small parameter L w /L = 
(uju L d )- 1 / d = (A/lu) 1 ^: 

P(uj, q) = {-iu)- 1 [F (z) + (A/-^) 1 " 7 , (19) 

where the scaling function $(z) has the same analytical properties as F(z). 

To find 7 one uses Eq. flI2p in the limit L u q <C 1. Substituting there D(u) oc 
I?~ d 1 + (L u /L) 1 / U (resulting from the standard renormalization group equation) instead 
of Eq. fll4l) , one expands the diffusion propagator up to the first power in D{ui)q 2 /ui. Com- 
paring such an expansion to Eq. (0), we have: 

7 = 1- (ud)-\ (20) 

Note that 1/2<7<1 due to the Harris criterion |J. 



Repeating the procedure which led to Eq. (|18D with P(u, q) given by Eq. (|T9|), we obtain: 



^(-l) n 5 d | X 2n + l| 2 A-^ f d% 



^2» = ™ 77 11 



/3vr 2 ci w " 7 (2tt 



x 



oo ' 

o C 2n IV— 



V J pikj cos 0j 



u 2n / j=l \ J 

Here, in contrast to Eq. fll8|), the integrand has an odd in z part. This is the only part 
which contributes to the integral (pip. As this integral is a nonzero dimensionless number, 
we obtain using Eq. (|20|): 

i?(s) = -c^/TV 2 ^ (s = w/A»l), (22) 

where Q/j is a numerical factor. For /3 = l,d = 2 + e expansion gives v — 1/e and 7 = 2/d 
near d = 2. In this case, the integrand in Eq. ( |2"T| ) has no odd part and q vanishes at d — 2. 

With .R(s) from Eq. fl22| ) the integral in the sum rule (||) is convergent, and we can 
use it for calculating the first integral in Eq. @. The second integral in Eq. (||) is also 
determined by the region s ~ (N) ^> 1, and we arrive at the announced result ([[]), where 
a d = 2^/7(1 -7). 

Since the coefficient must be positive, q > 0, and the correlator R(u) is negative for 
uj ^> A. For small cj^C A one can use the same zero-mode approximation Q as in the metal 
region for u <C 1/td, so that the correlation function R(s) should have the Wigner-Dyson 
form. We can conclude, therefore, that the energy levels are repelling at all energy scales. 

Note in conclusion that the Wigner-Dyson statistics can be represented as the Gibbs 
statistics of a classical one- dimensional gas of fictitious particle with the pairwise interaction 
V(s-s') = — ln|s — s'\. The Poisson statistics corresponds to V(s-s') = 0. If we suppose 
that the statistics of energy levels in the critical region can also be represented as a Gibbs 
statistics with some pairwise interaction V(s-s'), then such an effective interaction may 
be found, using the approach developed for the random matrix theory m^HTIf. Thus, in 
order to reproduce the asymptotics of the two- level correlation function (p2|), the interaction 



should have the form 115 : 
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V(s-s') = 




This interaction is valid for s — s' ^> 1. For small s — s' the interaction should be of the 
Wigner-Dyson form. Therefore, V(s — s') always remains repulsive. 

In order to check the conjecture about a pairwise nature of the effective interaction, one 
should investigate the higher order correlation functions. If they are factorizable like in the 
random matrix theory ||, then the Gibbs model with the interaction (p3|) will describe the 
whole statistics at the mobility edge. 
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